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ABSTRACT 

The  Coppens-Sanders  theory  for  the  one-dimensional,  nonlinear, 
acoustic  wave  equation  with  dissipative  term  describing  the  viscous 
and  thermal  energy  losses  encountered  in  a  rigid-walled,  closed  tube 
of  large  length-to-diameter  ratio  was  applied  to  finite-amplitude 
standing  waves  by  the  use  of  the  Fast  Fourier  Transform.  Computer 
programs  were  written  to  determine  the  amplitudes  and  phases  of  the 
first  255  harmonics.  Curves  of  harmonic  distortion  as  a  function  of 
the  strength  parameter  were  found  to  be  in  excellent  agreement  with 
available  experimental  data,  to  agree  with  the  Coppens-Sanders 
perturbation  analysis,  and  to  extend  the  theoretically  describable 
regime  closer  to  strengths  leading  to  the  formation  of  the  shock  front. 
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p,p  =  instantaneous  density,  ambient  density 

w    =  driving  (angular)  frequency 

u    =  resonance  frequency  for  the  fundamental 

r 

th 
c    =  phase  speed  in  the  tube  associated  with  the  n   harmonic 

1/2 
c~   =  ^3p/9p)      for  dissipationless  propagation  (linearized  theory) 

o  n— p 

M  Ho 
k    =  w/c 

L    =  length  of  the  tube 

x    =  Lagrangian  coordinate  measured  in  the  x  direction  from  the 
driven  end 

5  =  particle  displacement  in  the  x  direction 

u    =  particle  speed  in  the  x  direction,  (3C/3t) 

p    =  acoustic  pressure 

p   =  acoustic  pressure  for  the  n   harmonic 

P    =  amplitude  of  p 

2 
A    =  P  /Mp  c 
n     n   o  o 

<f>    =  phase  angle  associated  with  the  temporal  part  of  the  acoustic 
pressure 

<5    =  dissipation  parameter 

6  =  dissipation  parameter  for  the  n   frequency  component 

2 
UT   =  D ' Alambertian  in  Lagrangian  coordinates  (one  spatial  dimension) 

0)    =  dissipation  operator  for  the  n   frequency  component 


Ln 


(6  /u  )33/323t  _  6  92/9*2 
n  n    x     n 


M    =  a  Mach  number  (See  Eq.  2.3) 

2    2 

B/A  =  parameter  of  nonlinear ity  p  (3  p/3p  )  _ 

P  v. 


b         =1+1/2    (B/A) 

H         =  amplitude  operator    (See  Eg.   2.13) 


n 


=  phase  operator  (See  Eq.  2.14) 


2 

B  =  amplitude  associated  with  p  (See  Eq.  2.9) 

r  =  phase  angle  associated  with  B 

£  =  (poCo2)_1  (nB  +  4n/3) 

nD  =  coefficient  of  bulk  viscosity 

a 

n  =  coefficient  of  shear  viscosity 
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I.   INTRODUCTION 

Theoretical  investigations  of  finite-amplitude  waves  have  been 
characterized  by  various  approximation  and  perturbation  techniques, 
but  more  significantly  by  the  postulation  of  various  models  of  the 
absorptive  mechanisms [1,2,3,4] .  Coppens  and  Sanders  [5]  have  in- 
vestigated finite-amplitude  standing  waves  for  rigid-walled  tubes 
of  large  length-to-diameter  ratio  by  a  perturbation  method  and  have 
proposed  that  the  dominant  energy  losses  are  produced  by  the  viscous 
effects  of  the  fluid  at  the  walls  and  the  thermal  losses  due  to  the 
wall  material.  They  have  shown  that  an  approximate  wave  equation  for 

the  process  is 
oo 


DL+  °\n  KV 


where  UL  —  )  ^-v\    and  the  subscript  n  represents  the  n   harmonic 
of  the  waveform.  The  n   frequency  component  of  0.      is 

■HVi 

where  6  is  the  dissipation  parameter  [3,5]  for  the  n   frequency 
component  defined  by  6  =  6  / '    /- .  The  process  has  two  characteristic 
parameters,  one  related  to  the  strength  of  the  standing  wave  and  the 
other  related  to  its  frequency. 

The  purpose  of  this  thesis  is  extend  the  Coppens-Sanders  approach 
by  a  "Fast  Fourier"  analysis  technique.  By  an  iterative  method,  the 
amplitudes  and  phases  of  the  various  harmonics  are  obtained  for  dif- 
ferent values  of  the  strength  and  frequency  parameters  and  are  compared 
with  available  experimental  and  theoretical  results. 
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II .   THEORY 


In  this  section  a  Fourier-synthesis  approach  is  presented  for 
the  one-dimensional,  nonlinear,  dissipative,  acoustic  wave  equation 
based  on  the  theory  of  Coppens  and  Sanders  as  applied  to  finite- 
amplitude  standing  waves  in  a  rigid-walled,  closed  tube  of  large 
length-to-diameter  ratio. 

We  begin  with  Eq.  1.1  which  describes  the  process  in  question. 


Notice  that 

oo 


Po 


C 


o 


so  that  Eq.  1.1  may  be  written  as 


OO 


W-\ 


nNadtfv   b    * 


z. 


L. 


\_r\ 


pGcoz    ^)XZ 


[& 


(2.1) 


(2.2) 


Assume  that  p  may  be  written  as  a  Fourier  series 
oo 


<P 


n-\ 


(2.3) 


where  P  =Mp  c  A  . 
n    o  o  n 

Substitution  forDT   and  CD,,      into  the  left-hand  side  of  Eq.   2.2 
L  Ln  ^ 


results  in 
OO 

2. 


n=\ 


^' 


,2 


y       .   On    ^> \ 


^       Ci^t2        H,     ^^ 


a* 


^y< 


"tPn*  (2.4) 


Upon  substitution  for  6     and  u   ,  Exp.      2.4  becomes 

n 


n=\ 


n*  ^     Co2-  Bt2  T  n^  V  3t 


v\ 


(2.5) 
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2.  "7  2  "P 

Notice  that  ~^2  -  _(n  j^l  -^n     and^g)  -  L  J  ^ 

or       \       /  ax       \     ' 


so  that  Exp.  2.5  may  be  written  as 

CO 


->  hjy 


,2 


Si  \ 


S,     ^ 


n2 


£ 


2^Z 


3 


& 


r\ 


(2.6) 


"c0"     n2^    ^"t 

Define  u  =co  +  Aco  where  Aw  <  <  u  ,  and  for  the  fundamental  JR  =  ^V  —  _ 


n-  \ 

ie  oi  =w     +  Aw  where  Aw  <<   c 
r  r 

where  the  approximation  is  valid  for  weak  absorption.  The  approximation 


may  be  made  that 


2 


CJ 


C 


2^2 

o 


re 


c, 


+ 


c> 


c, 


ZAgj 


(2.7) 


Define  (c,/c  )  =1-6  and  neglect  second-order  terms;  Exp.  2.6  may 


then  be  written  as 

■».   .2.r  . 


n 


w 


Lrv 


ZAcj 

wrS, 


3 


^P 


n 


(2.8) 


n-\ 


The  right-hand  side  of  Eq.  2.2  contains  the  quantity  p  .  This 

involves  multiplication  of  two  infinite  series.  Terms  of  the  form 

cos  p(L-x)  sin  (qwt  +  $  )  for  p  ?  q  may  be  neglected  [5] ,  so  that  p^ 

si 

may  be  approximated  by 

.£>  =    H|SJ  ^T  Br,cJosn-^(L-x)s\n(nojt+rn|  .  (29) 


n=\ 


Then  we  have 


^f(#)  '     ^2     2 


V* 


>X 


—  (\^^)-A)         and  Eq.   2.2  may  be  written  as 


s. 


-tH  + 


2Acj 


I      S 


n=v 


_p^        S  <*v    n^-cu  S*t_ 


AnCosn&fL-X  s'Jncjt-K^b.iO) 


oo 


=  ^2^  cosni^L-x  sWi (nut  -f-rn  J  < 


n=\ 
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If  attention  is  limited  to  the  rigid  end  at  x  =  L,  Eq.  2.10  may  be 


considerably  simplified  to 
2ACJ 


Ln 


•/2 


+ 


\ 


S.^r- 


sm(acot  +<f£j  -  ^2co^(n^t^4>  A 


n    (2.11) 


This  becomes 


1  rv=\ 


)  An^(^"t+(/)0+^) 


(2.12) 


n-\ 


where 


_  Mb 


o& 


Z^/^  H^Bn^^(n  ^  +rn) 


rm 


Hn" 


and 


Lvn 


/2. 


ZAoj 


,2 


2 


n 


Vz 


Z 


(2.13) 


^="t 


-I 

an 


n 


/z 


Vz 


ZAcj 


(2.14) 


The  parameters  relating  to  the  strength  and  frequency  response  of  the 

standing  wave  are  Mb/ 6.  and  — —   ,  respectively. 

r  1 
An  iterative  solution  to  Eq.  1.1  may  be  obtained  by  successive 

application  of  Eqs.  2.12,  2.13,  and  2.14  to  an  initial  set  of  A  's 

and  <j>  's.  The  necessary  procedure  is  shown  as  a  flow  diagram  in  Fig. 

1.  A  discussion  of  the  Fast  Fourier  Transform  and  the  reasons  for  its 

use  herein  are  found  in  Appendix  A.  Notice  that  any  D.C.  terms  which 

appear  in  the  procedure  have  to  be  set  equal  to  zero  since  the  theory 

excludes  these  terms. 
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A  source  of  computational  difficulty  arises  when  the  frequency 

2Ato 

parameter     g      becomes  unity.  For  this  case  Eq.  2.13  reduces 
r  1 


to  H  =\/n/2,  a  result  which  causes  the  computer  solution  to  become 

unstable  because  of  the  strong  excitation  of  the  higher  harmonics. 

The  reason  for  this  difficulty  may  be  seen  as  follows:  the  resonance 

frequencies  for  the  classical  nondissipative  case  are  defined  by 

w  =  no)  ,  where  w  =  c  —  .  For  linearized  propagation  in  a  duct, 
n     c        c    o  L  ^  *^ 

resonance  frequencies  for  exciting  each  harmonic  would  be  given  by 

w   '   =  noj1   =  n=r  c   ,  where  c     =  c      (1-6  /2)  .     As  n  ■*■  °°,  we  have  the 
n  Xi     n  no  n 

i 

limits  c  ->  c  and  a/  =  to  .  For  the  case  of  large  n  then,  since 
n    o  c  n  r 

a)  =  w  +  A  co,  the  value  of  the  response  parameter  corresponding  to  the 

2Aw 
predicted  resonances  in  the  linearized  theory  is  given  by    g   =  1; 

r  1 
thus,  the  higher  harmonics  are  optimally  excited  and  problems  will 

arise  in  the  computer  program,  which  admits  only  the  first  255  harmonics, 

Notice  that  the  bulk  losses  have  been  neglected  up  to  this  point. 

If  these  terms  are  retained©   may  be  rewritten  [5]  as 


© 


Ln 


s 


n 


U) 


+  £ 


n 


a 


'h 


-S. 


V 


^U^z       "n   ^x- 


(2.15) 


In  most  experimental  configurations  the  value  of  e  is  usually  so  small 
that  it  may  be  neglected  when  compared  to  6/0..  Injection  of  an 

unrealistically  large  value  of  e  does  alleviate,  to  some  extent,  the 

2Au) 


instabilities  in  the  neighborhood  of 
and  2.14  become 


n 


2 


n 


/z 


+ 


2Acj 


+ 


u)  6, 
r  1 


=  1 


Equations  2.13 


n 


y-2     / 


1 

2 


(2.16) 
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""  a  -+    "'    J/n/z  -en 


+  -T7,  + 

(2.17) 


n'*  '   "V  S 


and  H  would  have  the  asymptotic  value  of  1/e  n  as  n  ■*■  °°.  The 
effectiveness  of  this  modification  will  be  discussed  in  more  detail 
in  Section  III. 
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Initialize  A',tf>  ,  —r—  ,   — r—  ,   Correction  Factor 
n  n   o,   w  6- 


Calculate  H  ,  6 
-a: — n_ 


Calculate  P  =  A_ 
m    ti 


sin 


511   yn, 


m  =  0  to  512 


R  =  P* 
m    m 


i 


Fast  Fourier  Analyze  P^  to  obtain  Bn/r'n 


1 


A'  =  A  + 


M  b 


H  B   -  A 


6-.  n  n n_ 


Correction  Factor 


b  *  =  4  +  ( r  -9  -tj))  Correction  Factor 
n    n    n    n   Yn 


Punch  cards  and  print  output  for  A  ,  <j> 


h 


Non-carrputer  step:     Have  A     and  <f>     reached  stable  values? 


Flow  Diagram  for  Solving  Eq.  1.1 


Figure  1 
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Ill .   RESULTS,  DISCUSSION,  AND  CONCLUSIONS 

The  computer  program  written  to  solve  Eq.  1.1  using  Eqs.  2.12, 
2.13,  and  2.14  gave  results  consistent  with  the  theoretical  pre- 
dictions of  Coppens  and  Sanders  [5]  and  Ruff  [10] .  An  additional 
perturbation  solution  that  was  performed  to  insure  that  the  program 
was  operating  properly  may  be  found  in  Appendix  B.  An  initially 
sinusoidal  pressure  waveform  when  used  with  a  correction  factor  of 
unity  gave  rapid  convergence  to  a  distorted  waveform  except  that  in 

the  neighborhood  of  response  parameter   the  solution 

o)r61 

became  unstable.  This  region  of  instability  gradually  increased  as 
the  strength  parameter  was  increased.  For  these  cases  the  modifications 
defined  by  Eqs.  2.16  and  2.17  were  used  with  e  =  0.001  and  allowed  a  few 
more  values  of  P  /P,  to  be  computed,  but  was  not  successful  in  gener- 
ating stable  results  for  the  higher  strengths.  The  excellent  agree- 
ment with  Coppens  and  Sanders  is  illustrated  in  Fig.  2  for  the  strength 
curves.  Figures  3,  4,  and  5  illustrate  both  the  area  where  the  solution 
deteriorated  and  also  the  frequency  characteristics  of  the  various 
harmonics. 

Since  the  only  difficulties  encountered  were  in  the  vicinity  of  a 
response  parameter  of  unity,  it  seems  plausible  that  the  instability  is 
a  result  of  the  behavior  of  the  higher  harmonics  which  became  important 
in  this  region.  Recall  that  a  response  parameter  of  unity  corresponds 
to  the  optimum  excitation  of  the  higher  harmonics  and  that  of  necessity 
a  truncated  series  has  been  used.  Then  the  higher  harmonics  that  are 
retained  are  affected  by  the  absence  of  the  harmonics  that  have  been 
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discarded  and  react  by  growing  in  amplitude  to  absorb  the  energy  that 
would  have  been  dissipated  by  those  terms  that  have  been  omitted.  The 
next  refinement  of  the  theoretical  development  of  this  research  would 
be  to  cause  the  retained  harmonics  to  behave  as  if  the  neglected  higher 
harmonics  were  present. 

Since  a  solution  using  iterative  processes  is  easily  adapted  to  a 
computer  simulation  it  is  suggested  that  this  method  be  continued.  The 
use  of  a  time-sharing  system  is  recommended  since  it  allows  close  control 
of  the  program  and  saves  a  great  amount  of  time. 

It  has  been  shown  that  a  Fourier-analysis  approach  to  the  Coppens- 
Sanders  theory  for  the  one-dimensional,  nonlinear  acoustic  wave  equation 
is  readily  adapted  to  computer  simulation  and  extends  the  region  of 
validity  to  the  pre-shock  regime.  The  series  truncation  appears  to  be 
the  only  obstacle  to  obtaining  a  solution  valid  for  the  entire  spectrum 
of  strength  and  response  parameters. 
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Figure  2 
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/     o  Without  bulk  absorption 
■  With  bulk  absorption 


Response  Curves,  Mb/6  =0.2 
Figure  3 
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Response  Curves,  Mb/6.  =0.4 
Figure  4 
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Response  Curves,  Mb/6,  =0.6 


Figure  5 
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o  Without  bulk  absorption 


Response  Curvw,  Mb/6.  =0.8 
Figure  6 
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APPENDIX  A 

The  Fast  Fourier  Transform 

The  Fast  Fourier  Transform  [6,7,8,9]  provides  an  effective  means 
for  estimating  the  Fourier  coefficients  of  a  set  of  evenly-spaced 
discrete  data.  The  Fast  Fourier  coefficients  differ  only  slightly 
from  the  conventional  Fourier  coefficients.  Construction  of  wave- 
forms is  of  prime  interest  for  the  higher  strength  parameters,  where 
the  Fast  Fourier  Transform  has  the  distinct  advantage  of  eliminating 
the  Gaussian  overshoot  characteristic  of  the  conventional  Fourier 
Transform.  This  property  is  illustrated  in  Figs.  7,  8,  and  9.  The 
waveform  of  Fig.  7  was  sampled  to  provide  512  evenly-spaced  discrete 
data.  This  was  the  input  to  the  procedures  available  on  the  IEM-360 
Computer  that  perform  Fourier  and  Fast  Fourier  analyses.  The  coef- 
ficients from  the  Fourier  analysis  were  used  to  construct  the  waveform 
of  Fig.  8,  while  those  of  the  Fast  Fourier  analysis  were  used  for  Fig. 
9.  From  the  standpoint  of  computer  utilization,  the  Fast  Fourier 
Transform  has  the  advantages  of  requiring  less  space  and  less  time 
and  was  used  for  the  analyses  in  this  research. 
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?C2 


Waveform  for  Analysis 


Figure  7 
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Reconstructed  Waveform  Using  Fourier  Transform 


Figure  8 
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Reconstructed  Waveform  Using  Fast  Fourier  Transform 


Figure  9 
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APPENDIX  B 


Perturbation  Solution  of  Eq.  1.1 

Assume  as  a  first  approximation  that  the  solution  to  Eq.  1.1 
consists  of  only  the  fundamental  which  may  be  written  as 


-(pi  =  p0C0  Mcos-teCL-^smu/t  . 


(B.l) 


Then 


6>M^ 


COS2Jr(\— x)  -COSZc^t 

-coszx(l-x)  cosecut 

so  that  the  right-hand  side  of  Eq.  2.2  is,  at  x  =  L, 


(B.2) 


b     ^IMl-  hp0cZMi2p+cos2ajtl.  *-3) 

n  r2   ~\  Z  -   r o  o       L  J 

Po^o    OX 

The  squaring  of  the  sin(ut)  gives  a  D.C.  term  and  cos(2wt) . 
The  D.C.  terms  may  be  excluded  [5]  so  that  the  second  harmonic  required 
on  the  left-hand  side  of  Eq.  2.2  may  be  assumed  to  be  of  the  form 

£pz=  Pz  cosZ$(L-x)s^(zojt-^z)  .  (b.4) 

Then  -(D    +£>  9  )-fp2.    mSiY  ^  written  as 

+  2VF  ^2S ,  cos2&  (v.-  ^)cos(2cJt+^<B-  5> 

For  Aw  =  0,  .LO  o     so  that  at  x  =  L,  Exp.  B.5  may  be  written  as 


-4?z9t8,  "I  +/vT  5'm(Zujt^)-|/-VT|co5(2a3t+4); 


(B.6) 
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which  when  simplified  is 
2 


Z- 


-4R  Jk  S,  .7G8s\n(zcjt+c/)o+247.S0 


-  o 


Equating  Exps.  B.3  and  B.7 

-4  Pz^S,  [.-7 (08 -5mfecjt+^4-2^-/„b 

2      2^Zr 


-bp0c0  Mi^  L"|4"eo<52-wt| 


(B.7) 


(B.8) 


Mb 
which  may  be  solved  for  P„  and  <j)_.     For  —  =  0.1, 

?z-  .0325p0c0ZM      and     <£z^  .39racUcms  . 

This  agrees  with  the  predicted  values  of  Coppens  and  Sanders  [5] 
and  Ruff  [10]  and  the  results  of  the  computer  program  written  to  solve 
Eq.  1.1. 
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COMPUTER  PROGRAM 

This  section  contains  the  computer  program  used  to  solve  Eq.  1.1 
and  a  list  of  the  symbology  pertinent  to  the  program. 


AAMP(N) 

= 

A  i 
n-1 

FEE(N) 

= 

Yn-1 

BAMP(N) 

= 

Bn-1 

GAMMA  (N) 

= 

rn-l 

H(N) 

= 

H  i 
n-1 

THETA(N) 

= 

0  i 
n-1 

P(N) 

= 

amplitud 

R(N) 

= 

amplitud 

ARG 

= 

wt/511 

RESPM 

= 

'  r  1 

STRPM 

= 

Mb/6 

CORR 

= 

Correcti 

amplitudes  of  p/Mp  c   for  increment  N 

2  2 

amplitudes  of  (p/Mp  c  )   for  increment  N 


Computer  Program  Symbology 
Table  I 
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Computer  Program  for  Solving  Eq.  1.1 


PIM^N^TPN    A?Wp(2K^)  ,  CFP  ^"M  ,  B^  M0  f  2^7  )  ,G~' MMM  257  )  t 
1H(2C7) ,THFT^ (2C7) f TMV(fA) ,^(54)  fP(60C)tB fft^O) 
PJ  =  3.  I4le<33 

apn=2*P!/5ll. 

r 

r 

f     tmtt!M!7F    pC«:om,ctppm    fnp.Rf  aamo(N)  fF«==f  N). 

C 

r 

PFJ.n    (5,2)    orSpM»5TRPWfcnpo 
2    FPPM^T(3P6.3) 

PCS.D(5,*)     (MMP(T  )  ,P^C  (I)  ,1=1  ,2C7> 
a    copm* T( 3(^xt2PlC.7) ) 

C 

r 

r    tHj<:    crriin^    C^L^UL^TF^    tmc    /•MDLTTiinc    rjpFRATrRS    H(N)    ANT 

r  tHc     ph£^f    nrppftTpsc     T^=T1(N), 

r 
r 

nn   l    1=2,256 

J  =  !-l 

XN=J 

perips-f i.o/SO»t(xn) ) 

ni)MMYl=-Pft  Tl  P-l  ,04-RFrow 
C 
C  cqii^tjpk    2.13. 

H( T )  =  1.C/SQRT(DUMMY1 *0UMMY1+P r TP*° AT  if ) 
C 
C  FOUATiPN    2.14. 

TMFT*  (  I  )=ATA  N2(P6Tin,nilMMYl  ) 

r 

1     rnN'TINt)F 

P°    10Q0    \id*«:c  =  i  f27 

r 

r 

C     THIS    *cf,TIPN    C4LCUIATFS    THC    AMO|_j  TUDE*  «     BAMP(N)     *NIO    THF 

C  nH*SP^    r,«MM4fM)     CPRPFcCPMni  NG    Tp    TH^    ^0U4PF    PF    TMF 

C  PDFcqjpc, 

r 

r 

DP    11     1=1, c12 

off  )=0.C 

nn    12    J=l,257 

°(  T  )=P(I)+AAMP<  J)*SIN(ARG*U-1)  *(I-1)+FEC(J)) 

12  CnNTIN1IF 
R(IJ=(P(T)*P(I)) 

11  CPNTTNUF 
C 

r  *  C    TN!D|iT.     Q(i)     fPNT^TN'^     THC    AMCLITUPfS  ,     P0&     INCRCMENT 
f  "T  ,  «    CrtJ0F$PPNPT  NG     TP    Tmc    ^qii^dc    PF    THF    DP<=Srt)PF. 

C 

FALL     PM/pm    (PfPfTNV,c,TFFRP) 
C 

r  &S    niiTDUT,    p(i)    CCNT.SINS    tmc    p4«;T    FniJRjEp    COEFFICIENTS 

C  HF    THF    <:QUA0P    PF    tmf    PPF$SijRe:     (A0)/2,     BC=1,     U,     Bi,.. 
C  (&N)/2,    BN=0. 

r  tmc    ANjtc    ipp    Thf    CDEFFICTPNTS    nF   THF    CnSTNc    TcrW<;    /^in 
C  THF    PN«<    APF    the    CFEFPICIFNTS     PF    THE    « INF    tfpmc, 

r 

nn    13    1=3,511,2 
f  WP  t  tc     c££H    rCMPPMFNT     PF     ^w  r    «;qh/\rc     pp    THF    orf^cijrf    /c 

C  RSmd     cjm(  M(  i  NiCPFMENT)        +        G£"M*(N))     BY    II^ING    THE    PArt 

r  fpuptfr  rnppcir i^mts . 

J  =  f  T  +  l) /2 

RAMP(J)=<QRT(R(I*l)*P(I*l)+RfI)*P(T)) 
GAMMA ( J)=ATAN2CR(I) ,R(I«-1) ) 

13  CONTINUF 
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c 

C  SFT    n.C.    TPRVS    TC    7FRC. 

C 

RAMD(  l)=C.O 

HAMMA  {  1)=0.0 
HAMM/i  (  ?c7)  =Q.O 

c 
c 

C    this    ^FCTTON    ADJUSTS    THE    LFFT-HANT    SIOF    OF    EOIMTION    2.12 

C  TQWAPP?    THF    RICHT-H*ND    SIOF. 

f 

r 

nn    22    1=3,256 

ASMD(T)=/iiMp(j)+cnRR*(RAVP(T)*<;TQOM*Hn)*0.5-AAMP(T)) 
Fcf(  j  )=FFF(I  )♦<  GAMMA  (I  )-PPF(  T)-TutcT&(i  )  )*CTRR 
22    C0NTTNUc 

r 

C  SFT    n.f.     TPPMS    TQ    7FOC. 

r 

A  A  Mr ( i)=c.n 
ppcI 1)=C.C 
/i  «MP(  2S"')=C.O 
FPc(?57)=C.O 

C 

r 

r  PP!NT    tNH    PUNCH    0UTP»Jx. 

C 

r 

WRTTF    (6,3A)    NPASS 

34  FnpM«T(«     PAS*    MIJMRFR     *  ,!F) 
nn    32    1=1,10 

J  =  T-1 

WPTTF(  ^,35)      J,A£MP(T  )  tF<=p(  J)      B.»MO(  I  )  fGAMMA(  I  ) 

35  cnPMAT    <!5,4F1R.7) 

32  cnMTTMiF 

On    31    1=11,257,10 

J  =  !-l 

WOTTF     (6,33)     J,AAMP(T  )  ,P«=F(!) 

33  FnpM;T(Tct2P1<;(7) 

31    CONTINUF 

lOOH    cnNTTKnir 

WOITF     (6,45)     RFSPM,STRP^tCnoR 
4  5    cnPM/  t(  io X,  »PF5PM=t  ,pa,3  ,5X,«STRCM=«  ,  P^  .  3,  SX  ,  •  r.PRR=  •  , 
1P5.3) 
tMJNCH    43     ,(I  ,AAMP(  T)  fFEEU  )  ,1=1  ,257) 
43    ppdmaT(3( I^,2F1C.7) ) 

CTHD 

fnd 
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